#include <stdio.h>
#include <math.h>


double *geary_test  (double *x, int n)
{
  int i;
  static double y[2];
  double diff, s=0.0, mean = 0.0;

  y[0] = 0.0;
  for (i = 0; i < n; ++i)
    mean += x[i];

  mean /= n;

  for (i = 0; i < n; ++i)
  {
    diff = x[i] - mean;
    y[0] += fabs (diff);
    s += diff * diff;
  }

  s *= n;
  y[0] /= sqrt (s);
  y[1] = (y[0] - 0.7979) * sqrt ((double) n) / 0.2123;

#ifdef NOISY
  fprintf (stdout,"  TEST2  GTN    =%10.4f   Z(GTN) =%10.4f\n", y[0], y[1]);
#endif				/* NOISY */

  return y;
}
